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MODELING  m  SIMULATION,  TESTING  m  VALIDATION 


•  Outline: 

>  An  idealized  concept  of  a  V-hull  vehicle  design  for 
blast  analysis  has  been  studied  in  two  different 
commercial  software  packages  StarCCM+  &  LS-DYNA 

StarCCM+  is  :  Eulerian  ;  Finite  Volume 
LS-DYNA  is:  Lagrangian  /ALE;  Finite  Element 
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MODELING  m  SIMULATION,  TESTING  m  VALIDATION 


•  Introduction: 

A  section  cut  from  a  generic  3-D  model  is  used  study  the 
pressure,  velocity  and  temperature  distribution  due  to  the 
blast  effect  under  the  vehicle 


Section  Cut 
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MODELING  m  SIMULATION,  TESTING  m  VALIDATION 


•  Model: 

A  2-D  vertical  cut  has  been  made  from  the  3-D  vehicle  and  the 
following  Initial  and  Boundary  Conditions  have  been  applied  to  the 
model 
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Model: 

-  StarCCM+ 

-  Initial  Conditions 
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MODELING  m  SIMULATION,  TESTING  m  VALIDATION 


•  Model: 

-  LS-DYNA  Initial  Conditions 

-  Initial  loading  condition  required  a  conversion  from  total  pressure  to  the 
relative  volume  using  the  ideal  gas  law 

VfStands  for  Relative  Volume)  =  Initial  density  /  density. 

And 

Density  =  Pressure  /  (Cp-Cv)  *  T  (stands  for  Temperature) 

Cp  =  Specific  heat  under  constant  pressure 
Cv  =  Specific  heat  under  constant  volume 

For  Air  during  1  ms  to  10  ms,  the  value  gets  computed  by 

a  Load  Curve  input : 

Density  =  100  Bar/ 286  *  800  Kelvin 
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Concept  V-  Hull  Floor  Configuration  in  Two  Dimensions 


•  Mesh: 

3  levels  of  mesh  were  created  for  both 
StarCCM+  and  LS-DYNA  models: 

>  45,000  element  (Cell) 

>  145,000  element  (Cell) 

>  1,145,000  element  (Cell) 
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•  Analysis: 

Explicit-Unsteady-Non-Turbulent  Analysis 
-  Both  StarCCM+  and  LS-DYNA 

Courant  number  of  about  0.5  on  both  cases 
At  3  mesh  levels 
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•  Analysis  Cont’: 

-  Line  segments  on  the  Floor  and  the  Roof  were  used  to  compute  the 
resulting  forces  for  comparison  between  cases 
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MODELING  m  SIMULATION,  TESTING  m  VALIDATION 


•  Results: 

-  Lift  forces  on  the  Floor  and  the  Roof  were  compared 

-  Qualitatively  the  results  are  comparable 

-  Quantitatively  the  peak  and  residual  forces  are  distinct 


Floor  Force 
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•  StarCCM +  Results: 


Absolute  Pressure  (bar) 

0043416  L4 126  2J6  19 _ 4.F5H _ 6.5203 _ 6  3395 


Pressure  Distribution  StarCCM+ 


Temperature  Distribution  StarCCM  + 


Velocity  Distribution  StarCCM+ 


Mach  No.  Distribution  StarCCM+ 
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•  Result  Comparison: 
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Star-CCM+  Pressure  Distribution  LS-Dyna  Pressure  Distribution 


Star-CCM+  Velocity  Fluctuations 


LS-Dyna  Velocity  Fluctuations 
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•  The  Effect  of  Courant  Number: 

-  Few  different  time  steps  have  been  tried  to  investigate  the  effect 
of  the  Courant  Number 
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Courant  Number  Effect:  LS-DYNA 
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•  Courant  Number  Effect:  StarCCM + 
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•  Courant  Number:  StarCCM+  vs  LS-DYNA 
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MODELING  m  SIMULATION,  TESTING  m  VALIDATION 


Initial  Condition  Effect: 

>  The  static  supersonic  pressure  was  reduced  from  12.7  Bar  to  1  Bar  in 

StarCCM+  to  match  the  Total  pressure  of  100  Bar  in  both  StarCCM+  and  LS- 
DYNA  as  input 


Supersonic  Static 
pressure  12.6  Bar 


Supersonic  Static 
pressure  1  Bar 


Static  Supersonic  Pressure  Pulse 
12.6  bar 


Static  Supersonic  Pressure  Pulse 
1  bar 
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•  Supersonic  Static  Pressure  : 

LS-DYNA  vs  StarCCM+ 

Floor  force  Vs.  Time  when  Supersonic  Static  Pressure  is  set  to  1.0  Bar 
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•  Mesh  Density  Study: 
•  StarCCM+ 
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Floor  Force  Comparisons 
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•  Mesh  Density  Study: 

•  LS-DYNA 
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•  Full  FSI  Analysis: 

Conducted  only  with  LS-DYNA 


Flow  interaction  with  the  structure 
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Ful[  FSI  Analysis : 

LS-DYNA  computation 
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Possible  Future  Direction : 

>  Investigate  Blast  Analysis  in  different  Media  (Water/  soil) 

>  Compare  Other  FSI  software  against  the  current  one 

>  Import  Pressure  from  STARCCM+  as  a  1-way  Couple  to  LS-DYNA 
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Conclusion: 


Several  possibilities  have  been  investigated  that  include  studies  with  Courant 
number  setting,  Supersonic  Static  Pressure  setting  that  shows  Finite  Volume 
and  Finite  Element  results  can  be  equivalent 

The  Effect  of  Mesh  Density  on  Force  prediction  has  been  demonstrated 

A  2-D  type  of  analysis  has  been  demonstrated  to  mimic  that  of  a  full  3-D  analysis 
and  are  effective  for  what-if  scenarios 
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